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1. THE FORMULAS. In [3], [4] the first author gave a formula for the 
«th prime number that involves only the elementary operations 
+,-,x,4- and the floor function [ J. His conditional proof assumed certain 

inequalities based on the Prime Number Theorem. In this note, we use 
inequalities of Rosser and Schoenfeld [1] to give a complete proof of the 
slightly modified formula 
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, n>l. 



where 7i:(x ) , defined as the number of primes not exceeding x, is given by 
the formula 
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After the proof, we indicate various ways to modify and implement the 
formulas so that they operate in times 0((/zlog/z)^^^) and 0(^^^), 
respectively. 
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2. PROOF. For n a positive integer, let 
(3) 



J(n):=^l 

d\n 



denote the number of divisors of n. In [2], the first author found the 
formula 
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which holds since the quantity in parentheses is 1 or according as i does 
or does not divide n . 

Let F be the characteristic function of the set of prime numbers 



(5) 



F{n) := 



1 if n is prime, 
otherwise. 



From (3), we have d{n) = 2 if « is prime, and d{n) > 2 if « is composite. 
Since 2 < d{n)< n for n > 1 , we have the formula 



(6) 



F{n) = 1 + 



2 - d(n) 



n 



n>\ 



(7) 



Using (5), we write the prime-counting function k{x) as the sum 

J-2 



with the convention that any sum J^^-a ■^^^^ a>b. From (7), (6), (4), 
we obtain formula (2) for k{x ) . 

In order to derive formula (1) for from (2), we require the 
following lemma. 
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Lemma 1. For n>\ ,we have the inequalities 



(8) 



7r(2n\ogn + 2) <2n. 



(9) 



<2n\ogn + 2. 



Proof. Rosser and Schoenfeld [1] proved that 



(10) 



>n\ogn, all n. 



(11) 



<n\ogn + n{\og\ogn-l/2) , n>20. 



From (10), we have >2nlog2n. Since 7t(p2n) - 2n, it follows that 
K{2n\og2n) < 2n, which implies (8) if « > 1 . 

To prove (9) for n > 1 , we verify it numerically for « = 2,3,...,20, 
and note that (11) implies (9) for « > 20. This proves the lemma. • 

For n>\ , Lemma 1 implies that 



The desired formula for the nth prime number follows immediately. This 
completes the proof of (1) and (2). • 
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3. IMPROVEMENTS. As they stand, the formulas for 7r(x) and 
operate in times 0((n\ognf), O(x^) and 0(n), respectively. We can 
improve these bounds by modifying the formulas as follows. If i divides n , 
so does n/i, so to compute it suffices to consider only i<4n. This 
reduces the times for d{n) and k{x) to O^n^^) and 0(^/^), respectively. 
For p^, compute d{j) (and thus F^j")) for j<2n\ogn+2, then compute 
7t{k) recursively as 7t{k) = 7l{^-l) + F{k). This reduces the time for 
to 0{{n\ognf^^). 

We can also improve the computation times (but not the 0(.) 
bounds) in the following two ways. First, instead of the floor of «//, use the 
integer quotient of n by Second (as P. Sebah [5] has pointed out), in 
formulas (2) and (7) for 7Z{x), after j=l we only need to sum over odd 
numbers, after y= 3 only over numbers relatively prime to 6, and similarly 
for other moduli 30, 210, 2310,. . .,m. This "sieving" cuts computation 
time by a factor of m/(p(m), where (p is Euler's totient function. 

We thank P. Sebah for discussions on the bounds, C. Rivera for the 
square root optimization, and J. McCranie for the quotient acceleration. 
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